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I.  SUMMARY  OF  PROGRESS 

During  the  first  year  of  the  contract,  theoretical  methods  of  predicting 
ordering  parameters  of  carbon  metal  based  systems  have  been  Investigated. 
Preliminary  methods  of  calculating  the  critical  ordering  temperature,  and  the 
maximum  degree  of  order  have  been  used  to  examine  the  characteristics  of  C-Zr, 
C>Mo,  C-T1  and  C-V  systems.  Based  upon  these  calculations,  the  C-T1  system 
has  been  chosen  as  the  most  promising  system  In  which  the  ultrahigh  strength, 
ductility  and  temperature  resistance,  properties  desired  for  space  power 
generation  materials,  can  be  obtained. 

Fifteen  titanium  based  alloys  have  been  manufactured  using  an  arc  melting 
furnace  to  compare  experimentally  the  effect  of  carbon  content  on  the  ordering 
parameters  of  carbon  metal  alloys  with  theoretical  predictions.  The  alloys 
produced  contained  between  0  and  53  atom  percent  (0  to  22  weight  percent)  car¬ 
bon.  This  first  group  of  alloys  did  not  contain  transition  and  rare  earth 
additions  required  for  ductility  Improvement.  This  will  allow  an  unperturbed 
comparison  of  ordering  parameter  theory  with  experimental  results.  Optical 
metallography  and  hardness  tests  have  been  completed  on  these  alloys  with 
Interesting  results  for  alloys  containing  betwen  35  and  53  atomic  percent  (12 
to  22  weight  percent)  carbon.  These  alloys  were  found  to  contain  significant 
amounts  of  phases  not  predicted  from  the  phase  diagrams.  X-ray  diffraction 
tests  are  being  conducted  to  Identify  these  constituents. 

An  additional  ten  alloys  have  been  manufactured  using  a  new  Induction 
generator  purchased  under  the  contract.  Two  compression  tests  have  been  con¬ 
ducted  on  specimens  cut  from  these  melts.  One  of  these  specimens  was  a  T1-5 
w/o  carbon  alloy  which  was  found  to  have  a  yield  strength  of  100  ksl,  ultimate 
tensile  strength  of  164  ksl  and  5.4%  strain  to  fracture  at  room  temperature. 
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This  alloy  has  a  melting  temperature  of  nearly  2000**C  compared  to  traditional 
Ti  alloys  such  as  T1-6A1-4V  which  has  a  melting  temperature  of  1650*C,  UTS  of 
130  ksi  and  elongation  of  4%  in  the  forged  bar  condition. 


The  project  is  currently  on  schedule  and  within  budget. 

II.  OBJECTIVE 

The  objective  of  the  proposed  research  is  to  investigate  a  new  class  of 
materials  composed  of  30-60  atomic  percent  carbon  (C-Mo,  C-Zr,  C-V  C-Ti)  for 
ultrahigh  temperature  applications  in  space  power  systems.  The  proposed  alloy 
systems  exhibit  melting  temperatures  in  excess  of  2500**C  and  form  long  range 
ordered  structures  which  are  expected  to  yield  materials  with  exceptional  high 
temperature  strength  and  irradiation  resistance. 


III.  INTRODUCTION 

The  generation  of  power  In  space  by  methods  other  than  direct  conversion 
(photovoltaic)  will  require  new  concepts  In  materials  to  obtain  optimum  per¬ 
formance.  Traditional  materials  used  in  the  construction  of  land  based  power 
systems  such  as  zircaloy,  stainless  steel,  and  In  some  exotic  applications,  Nb 
and  V  steels,  do  not  have  the  necessary  combined  ultrahigh  temperature 
strength.  Irradiation  resistance,  and  strength  to  weight  ratios  that  are 
required.  The  space  reactor  will  operate  In  a  near  perfect  vacuum,  elimi¬ 
nating  oxidation  problems  which  plague  many  high  temperature  materials  (Mo,  V, 
etc.).  However,  sublimation  of  the  materials  (vaporization  at  high  tem¬ 
peratures)  will  become  a  controlling  factor.  In  space,  the  rejection  of  heat 
must  occur  by  radiation  of  energy  to  the  environment,  requiring  temperatures 
of  1000*C  or  greater  for  thermodynamic  efficiency.  Additionally,  the 
materials  are  required  to  withstand  the  radiation  of  the  nuclear  core  (be 
resistant  to  Irradiation  Induced  densificatlon,  swelling  and  creep)  and  be 
lightweight  for  transportation  Into  space. 

Viewing  the  above  general  requirements,  carbon  based  materials  (C 
materials)  are  a  superb  candidate  for  development.  C  fiber  -  C  matrix 
materials  have  been  highly  successful  on  the  space  shuttle  where  they  have 
been  used  for  components  ranging  from  cargo  bay  doors  to  high  temperature 
motor  casings.  Such  materials,  however,  are  not  suited  for  applications  which 
require  exposure  to  nuclear  radiation.  Testing  In  the  high  temperature  gas 
cooled  reactor  program  has  Identified  two  problems  with  C  materials. 

The  first  Is  the  permeability  of  the  C  material  to  radioactive  fission 
gases  [1,2].  Despite  the  development  of  ultrahigh  density  C  materials,  spe¬ 
cial  C  coatings,  and  differential  pressurization  of  components,  leakage  of  C 
clad  fuel  occurs  after  very  short  neutron  exposures  ((^t  <  1  x  lO^^n/cm^). 


Heat  pipe  components,  fabricated  from  pure  C  or  C  composite  materials,  would 
therefore  be  unable  to  maintain  proper  partial  pressures  for  operation.  The 
second  problem  results  from  volumetric  changes  (densificatlon)  that  occur  In  C 
materials  upon  neutron  Irradiation  [1].  These  changes  are  highly  orientation 
dependent.  The  C  fibers,  which  contain  material  with  the  basal  planes  aligned 
perpendicular  to  the  stress  axes,  undergo  a  substantially  different  volumetric 
change  than  the  matrix  material,  leading  to  the  generation  of  high  Internal 
stress,  cracking  and  ultimate  failure.  A  class  of  materials  which  promise  to 
overcome  the  above  described  problems  Is  ordered  C-metal  alloys.  Ordered 
alloys  offer  many  potential  advantages  over  conventional  alloys  at  elevated 
temperatures  [2-5].  The  atomic  ordering  produces  a  pronounced  Increase  In  the 
work  hardening  [5-9],  Improves  the  fatigue  resistance  [10],  and  retards, 
because  of  stronger  binding  and  closer  packing  of  atoms,  most  thermally  acti¬ 
vated  processes  such  as  creep  and  grain  growth  [11].  In  addition,  the 
strength  of  ordered  alloys  Is  less  sensitive  to  temperature  than  conventional 
disordered  alloys.  In  fact  some  alloys  show  an  Increase  rather  than  a 
decrease  In  strength  with  Increasing  temperature  up  to  the  critical  ordering 
point  [12].  The  critical  ordering  point  is  the  temperature  at  which  the 
material  reverts  from  a  defined  periodic  arrangement  of  alloying  elements  to 
the  random  arrangement  found  In  conventional  materials.  Electron,  Ion  bom¬ 
bardment  and  recent  neutron  Irradiation  results  show  the  long  range  ordered 
(LRO)  alloys  to  be  highly  swelling  resistant  [13-15]. 

Despite  the  above  advantages,  the  LRO  alloys  have  seen  only  limited  appli¬ 
cation  because  of  a  lack  of  ductility  associated  with  the  ordered  state 
[16-18].  Recently  Liu  and  his  coworkers  at  ORNL  have  succeeded  in  producing 
(Fe,  N1)3V  long  range  ordered  alloys  which  show  tensile  elongations  greater 


than  30%  [19-21].  This  has  been  achieved  by  controlling  the  ordered  lattice 
structure  through  use  of  the  e/a  ratio  or  average  electron  density  per  atom 
outside  the  Inert  gas  shell. 

Ordered  C-metal  alloys  represent  a  new  approach  which  combines  the 
desirable  qualities  of  C  (high  temperature  resistance,  low  weight,  low  vapor 
pressure)  with  the  radiation  and  permeability  resistance  of  ordered  materials 
These  alloys,  containing  30-60  atomic  percent  carbon,  all  exhibit  melting  tern 
peratures  In  excess  of  2500®C  [22,23].  By  careful  control  of  the  alloy  com¬ 
position,  desirable  ordered  structures  can  be  maintained  at  temperatures  In 
excess  of  0.7Tm  (approx. 1800*C) .  The  alloys  will  be  relatively  lightweight, 
yet  resistant  to  Irradiation  damage.  Control  of  the  ordered  structure  by 
modification  of  the  e/a  ratio  will  produce  compositions  with  good  ductility, 
strength,  creep  and  fatigue  properties.  In  addition,  the  C-metal s  will  ellml 
nate  the  permeability  problem  associated  with  pure  C  materials. 

The  objective  of  this  research  Is  to  Investigate  the  property-structure 
relationships  of  ordered  C-metal  alloys.  These  studies  will  provide  a  theore 
tical  and  experimental  basis  from  which  a  future  development  program,  leading 
to  optimization  and  detailed  characterization  of  ordered  C-metal  alloys,  can 
be  Initiated. 
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IV.  PROGRESS  DURING  FIRST  YEAR  OF  THE  PROJECT 

The  principal  accomplishments  during  the  first  year  of  the  project  deal 
with  both  theoretical  and  experimental  achievements.  Based  upon  the  initial 
proposed  project  plan,  the  research  is  on  schedule.  The  work  herein  reported 
was  performed  by  two  graduate  students  seeking  Master's  of  Science  and  Doctor 
of  Philosophy  degrees  from  our  Department  of  Mechanical  Engineering.  Most  of 
the  work  reported  was  performed  by  Pamela  Gills  who  began  research  on  the  pro 
Ject  in  June  1983.  S.  C.  Su  Joined  the  research  team  in  March  1984  and  as 
such  his  contributions  to  the  reported  progress  to  date  are  minimal. 

Approximately  twenty-five  percent  of  the  first  year's  effort  was  spent 
setting  up  theoretical  models  and  analyzing  the  C-Ti,  C-V,  C-Mo  and  C-Zr 
systems  for  the  alloys  with  the  highest  probability  for  success.  Fifty  per¬ 
cent  of  the  expended  effort  was  required  to  melt  and  examine  the  initial  C-Ti 
alloys  and  the  remaining  twenty-five  percent  was  used  setting  up  equipment. 
The  following  two  sections  describe  the  progress  made  to  date  in  the  theoreti 
cal  and  experimental  areas. 


i.  Introduction 


There  are  a  number  of  theoretical  models  available  which  allow  the  calcu¬ 
lations  of  the  critical  ordering  temperature,  maximum  degree  of  order  and  spe¬ 
cific  heat  as  a  function  of  temperature  for  ordered  alloy  systems.  It  is  the 
objective  of  this  research  to  apply  these  theories  to  the  Ti-C  alloy  system. 
The  results  of  the  above  computations  will  provide  a  theoretical  basis  upon 
which  the  most  promising  alloy  systems  and  ranges  of  compositions  can  be  cho¬ 
sen  for  melting  and  experimentation  under  /Uitivity  II.  It  is  hoped  that 
before  alloying  additions  are  made  that  we  will  be  able  to  extend  these 
theoretical  calculations  to  ternary  systems  [24]. 

This  section  documents  the  work  accomplished  so  far  in  these  efforts  and 
the  work  planned  for  the  future.  Included  in  the  report  is  a  brief  synopsis 
of  some  of  the  theoretical  models  of  ordering  found  in  the  literature,  and  a 
discussion  of  problems  encountered  in  applying  them  to  our  system. 

Unlike  conventional  alloys,  which  form  a  random  mixture  of  atoms  on  lat¬ 
tice  sites,  long  range  ordered  (LRO)  alloys  exhibit  periodic  arrays  which  form 
an  ordered  crystal  structure.  In  general,  unlike  atoms  tend  to  be  closer 
together  and  like  atoms  further  apart.  ITie  tendency  to  order  increases  with 
decreasing  temperatures,  therefore,  the  ordered  arrangment  of  atoms  must  be 
more  stable  or  energetically  favorable  at  lower  temperatures  than  a  random 
distribution.  For  this  to  be  true,  the  attraction  between  unlike  atoms  must 
be  stronger  than  the  attraction  between  like  atoms.  As  the  temperatures  of 
the  alloy  increases  the  energy  of  thermal  vibration  increases  and  when  this 
exceeds  the  so  called  ordering  energy  the  atoms  have  sufficient  energy  to 
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redistribute  themselves  Into  a  disordered  or  random  mixture;  it  Is  at  this 
point  that  the  disordered  structure  becomes  the  stable  configuration.  Thus, 
we  expect  a  transition  temperature  at  which  the  ordered  phase  becomes 
unstable.  This  transition  temperature,  the  critical  temperature,  T^,  is  ana¬ 
logous  to  the  Curie  temperature  in  ferromagnetism.  One  of  our  objectives  was 
the  prediction  of  the  critical  ordering  temperature  for  various  carbon  alloys. 
To  make  this  prediction,  a  review  of  the  various  theoretical  models  of 
ordering  was  undertaken.  This  review  is  summarized  in  what  follows.  Above 
the  critical  temperature  the  alloy  is  completely  disordered.  Once  the  criti¬ 
cal  temperature  is  reached  atomic  ordering  begins  and  at  each  temperature 
below  the  critical  temperature  (after  suitable  heat  treatment)  an  equilibrium 
degree  of  order  will  be  reached.  Thus,  at  the  stoichiometric  composition  of 
the  alloy,  the  equilibrium  degree  of  order  at  high  temperatures  will  be  zero 
and  at  0**K  the  alloy  will  be  completely  ordered.  Alloys  which  are  not 
stoichiometric  will  not  reach  perfect  order,  but  will  have  a  characteristic 
maximum  degree  of  order  depending  on  the  composition. 

The  first  suggestion  that  some  alloys  can  exhibit  an  ordered  structure,  or 
super-lattice  was  put  forth  in  1919  by  Tanvnann  [25]  after  his  observations  in 
the  Cu-Au  system.  Since  that  time  a  number  of  efforts  have  been  made  to 
theoretically  describe  and  predict  the  cooperative  phenomena  of  ordering  [26]. 

The  first  successful  theory  of  the  stability  of  superlattices,  as  a  func¬ 
tion  of  temperature,  was  developed  by  Bragg  and  Williams  [27].  Later,  Bethe 
[28]  developed  a  more  refined  method,  which  was  extended  by  several  authors. 
Neither  one  of  these  methods  were  very  clear  about  the  assumptions  which  were 
made;  Kirkwood  [29]  remedied  this  by  developing  a  theory  based  entirely  on 
the  standard  concepts  of  statistical  mechanics.  The  quasichemical  method  was 
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developed  also  at  this  time  by  Guggenheim  [30]  and  was  later  shown  to  be 
equivalent  to  Bethe's  method  [31]. 

The  more  recent  past  has  seen  the  development  of  several  more  theories, 
including:  the  Cluster>variation  method  developed  by  Takagi  [32],  Kikuchi 
[33],  Yang,  Li  [34]  and  Hill  [35],  the  method  of  concentration  waves  [11], 
and  a  band  theory  model  [37].  In  addition  to  these,  numerical  methods  based 
on  the  random  walk  or  Monte  Carlo  Method  [38]  have  been  developed. 

Except  for  band  theory,  all  of  the  above  models  require  the  interaction  or 
ordering  energies  of  the  alloy. 


(1) 

where  the  superscript  i=l  denotes  nearest  neighbors,  i»2  denotes  next  nearest 
neighbors,  etc.,  and  the  subscripts  denote  which  pair  of  atoms  the  energy 
represents.  Unfortunately,  this  energy  parameter  is  hard  to  come  by.  This 
will  be  discussed  later. 


11.  Bragg-W1 1 1 1 ams ,  Kirkwood,  Quasi  chemical  Methods 

In  order  to  specify  the  degree  of  order  of  the  atomic  arrangement  over  the 
lattice  sites  Bragg  and  Williams  Introduced  the  long  range  order  parameter  s 
In  their  original  paper  [27].  When  perfect  order  Is  achieved  In  a  binary 
alloy  of  components  A  and  B  the  lattice  sites  occupied  by  A  atoms  are  called 
a  sites  and  those  occupied  by  B  atoms  are  called  e  sites.  Let  N  be  the  total 
number  of  atoms  (A  and  B)  and  thus  the  total  number  of  lattice  sites,  (I.e.  no 
vacancies);  and  Fa  and  Fg  (where  Fg  »  1-Fa»)  be  the  fraction  of  A  and  B  atoms 
In  the  alloy  respectively.  Also  let  the  fraction  of  a  sites  occupied  by  A 
atoms  be  denoted  by  r^,,  (right  atoms),  the  fraction  of  o  sites  occupied  by  B 
atoms  be  w<j,  where  w^  *  l-r^  and  similarly  rg  and  wg  *  1  -  rg  represent  the 
rightly  and  wrongly  occupied  g  sites.  The  long  range  order  parameter  s  may 
now  be  defined  as 

f'a  ’  ^A  '*6  -  Pb 

s  - - - -  ,  (2) 

1  -  Fa  1  -  Fg 

Thus,  s  Is  defined  as  unity  for  perfect  order  and  vanishes  In  the  completely 
disordered  state. 

.  Now,  let  Q  be  the  total  number  of  all  types  of  pairs  In  the  lattice  and 
each  atom  be  surrounded  by  z  nearest  neighbors,  thus 

Q  »  (-)  N  .  (3) 

2 

Also  let  the  number  of  pairs,  which  are  AA,  BB  and  AB  be  denoted  by  Qaai  Qbb* 
QAg  respectively.  Consider  an  alloy  with  FaN  a  sites  and  FgN  s  sites,  then 
according  to  equation  (2)  the  average  nwnber  of  A  and  B  atoms  on  a  and  g  sites 
are  as  follows  [26]: 
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A  atoms  on  a  sites: 


rJFfiJi  =  (Fa  +  Fbs)FaN, 


A  atoms  on  b  sites:  wgFgN  »  (1  -  s)FaFbN, 

B  atoms  on  a  sites:  w^FaN  »  (1  -  s)FbFaN, 

B  atoms  on  8  sites:  TqFb  N  *  (Fb  +  Fas)FbN. 

The  number  g(s),  of  distinguishable  arrangements  or  configurations  for  a 
given  state  s  Is  then 

g(s) 


(4) 


CFaN]! 


C(Fa  +  Fbs)FaN]!C1-s)FaFbN]! 
CFbN]' 


(5) 


C(1-s)FaFbN]!EFb+Fas)FbN]! 

Assuming,  approximately  that  the  lattice  vibrations  are  Independent  of  the 
configuration  of  the  atoms  In  the  lattice,  the  partition  function  In  the  cano¬ 
nical  ensemble  may  be  written  as 


Z(T,s)  «  z  e-Wt/I^T  ,  (6) 

T 

where  T  Is  the  absolute  temperature,  s  Is  the  long  range  order,  K  Is  the 
Boltzmann  constant  and  Is  the  configurational  energy  of  the  state  t.  The 
summation  extends  over  all  the  states  t  of  a  given  order  s. 

Our  knowledge  of  Is  scant  at  best  at  this  time  and  Is  approximated  here 
as  consisting  mainly  of  the  Interaction  energies  of  pairs  of  nearest  neigh¬ 
bors.  These  energies  as  Identified  In  eq.(l)  are  Vaa»  Vbb>  Vab> 
corresponding  to  nearest  neighbor  pairs  Qaa»  Qbb»  Qab  respectively.  Therefore 
the  configurational  energy  Is  roughly  given  by 

Wt  ■  -  (Qaa  Vaa  QbbVbb  +  Qab  vab).  (7) 

for  any  specified  configuration  t. 

The  configurational  free  energy  F(s)  Is  given  by 


F(s)  »  -  kTiin(Z{T,s))  , 


(8) 


_ I.JmllUW  m  ■  lllllMtl»'l»IIW~  ■IIMIIIi  IMI  ■IIHIT- fc T  r  ^  »  r.  »  ^  ...«v 

^  • 

and  once  this  is  minimized  the  equilibrium  value  of  s  will  be  determined.  The 
corresponding  energy  of  the  alloy  is 

z  Wt  e-WtS 

E(s)  -  W{s)  «  -  ■  ^  -  .  (9) 

z  e“"TP 

T 

where  B  is  1/KT. 

Now,  using  the  conventional  procedures  of  statistical  mechanics,  we  define 
U(s}  by  the  relation 

^  e-»*TB  •  9(S)e-'J(s)B  (10) 

here  g(s)  is  the  same  as  in  eq.  (5).  Now,  8  can  be  written  as 

F(s)  »  U(s)  -  KTJin(g(s)),  (11) 

and  the  energy  of  the  alloy  can  be  written  as 


a  U(s) 

E(s)  - - 

3(1/T)  T 


(12) 


The  above  equations  for  the  free  energy  and  internal  energy  define  the  ther> 
modynamics  of  the  system  [26].  All  that  is  left  to  do  is  to  calculate  the 
partition  function.  This  however  has  proved  illusive  and  various  approxima¬ 
tions  have  been  made  to  it,  the  first  of  which  was  the  Bragg-Wi 1 1 i ams  approxi¬ 
mation  [27].  Their  approximation  consists  of  replacing  the  true  value  of  w^ 
for  each  configuration  of  s  by  the  simple  average  <Wt>  of  all  states  belonging 
to  a  given  s, 

<Wt>  -  -  {<Qaa>vaa  +  <QBB>  Vbb  +  <QAB>  Vab}.  (13) 

where  <QAA>t  <Qbb>  *'’‘1  <QaB^  t^*®  simple  averages  of  the  number  of  pairs  of 
each  kind  in  each  state  t  of  s.  Consider  <Qaa>»  AA  pairs  may  result  from  two 
arrangements:  (1)  one  A  atom  on  an  a  site  and  the  other  one  on  a  B  site,  and 
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(2)  both  atoms  on  s  sites.  On  the  average' there  are  FaN(Fa  +  Fgs)  a  sites 
occupied  by  A  atoms  each  having  z(1>s)Fa  neighboring  s  sites  occupied  by  A 
atoms.  There  are  then  zNFa^(1-s)(Fa  +  Fgs)  AA  pairs  of  type  (1).  The 
number  of  AA  pairs  with  both  atoms  on  0  sites  Is  ^zNFa^(Fb-Fa)(1-s) ^  [26]. 
Therefore, 

<Qaa>  -  2NFa^(1-s)(Fa  ♦  Fbs)  +  V2zNFa2(Fb-Fa)(1-s)2 
-  -  Fa^s^). 

and  similarly,  (14) 

<Qbb>  •  V2zN(Fb^-Fa^s^) 

<QAB>  ■  V2zN(2FaFb  +  2Fas^). 

Substituting  equations  (14)  Into  (13)  leads  to 

Nz  ,  .  , 

<Wt>  -  —  (fa^aa  +  2FaFbVab  +  Fb  Vbb) 

2 

Nz  ,  ,  (15) 

-  —  Fa^2Vab-VaA-Vbb)s% 

2 

Since  <Wj>  Is  Independent  of  temperature  In  this  approximation,  by  equations 
(10)  and  (12) 

<U^>  »  E(s)  »  U(s) 

Thus,  the  configurational  free  energy  becomes 

F(s)-F(0)  «  NKT[Fa(Fa  +  FBS)tn{FA(FA  +  Fgs) } 

+  Fb(Fb  +  F^s)tn{Fg(F3  +  F^s)} 

+  2F^Fg(l-s)&n{F^Fg(l-s)} 

-  2FA£nFA  -  2FB£nFB]  -  NzVFa^^ 
after  using  the  expression  for  g(s)  In  eq.  (5)  and  eq.  (15)  for  the  energy  In 
eq(ll),  where  V  Is  as  defined  In  equation  1,  In  the  nearest  neighbor 
approximation.  For  the  alloy  of  composition  Fa  ■  Fg  *  Vz 


(16) 


(17) 
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this  expression  reduces  to 


N 

F(S)  -  F(l)  -  -  KT{(l+s)tn(l+s)  +  (1-s) 

2 

tn(l-s)  -  2tn2}  +  V4NzV(l-s^).  (18) 

The  equilibrium  state  of  the  configuration  In  the  alloy  Is  found  by  minimizing 
the  free  energy.  Thus,  the  condition 


if 

—  •  0 

as 

zV 

gives  s  ■  tanh  - s  . 

2KT 

This  equation  determines  the  equilibrium  value  of  s  as  a  function  of  T. 
above  equation  can  be  decomposed  Into 


(19) 

(20) 


The 


2KT 

s  *  tanh  X  and  s  *  -  ,  (21) 

zNX 

When  these  tvfo  equations  are  plotted  as  a  function  of  X  the  Intersection  of 
the  curves  give  the  equilibrium  value  of  s  at  T  (see  Figure  1).  There  are  two 
solutions  of  (21),  one  at  s-O  and  one  at  s>0  for  all  T  less  than  some  Tc. 

The  solution  s»0  corresponds  to  the  maximum  free  energy  or  the  unstable  state 
of  order,  whereas  the  solution  of  s>0  corresponds  to  the  minimum  free  energy 
of  equilibrium  value  of  order.  The  solution  of  s>0  Indicates  the  existence  of 
the  ordered  state  at  temperatures  below  a  critical  value.  The  value  of  s 
decreases  very  slowly  at  first,  then  as  the  temperature  continues  to  Increase 
drops  off  abruptly  at  Tc,  and  vanishes  at  all  temperatures  greater  than  Tc. 

The  critical  temperature  Is  found  by  the  conditions 

3F  3^F 


3s 


s 


0. 
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(22) 


Figure  1:  Determination  of  S  at  T  from  equation  21. 


These  conditions  lead  to 


zV 

KTc  -  —  (23 

2 

for  the  critical  temperature  of  an  alloy  with  composition  Fa  *  Fb  *  Vz* 

A  higher  approximation,  developed  by  Kirkwood  [29],  Is  one  which  takes 
Into  account,  on  the  average,  the  energy  spread  of  for  a  given  s  around 
Its  average  value  <Ut>.  Although  It  Is  not  possible  to  know  the  actual 
distribution  of  we  can  find  the  moments  of  deviations  from  <Wt>.  We  write 
equation  (10)  as 


U(s)/KT  »  tn[g(s)]  -  in  [z 

T 


and  expanding  the  exponentials, 
U(s)/KT  »  -  ml  - 


<W^>  1  <W^S  1  <Wt' 


KT  2!  (KT)^  31  (KT)^ 


♦  mmm^ 


.  g(s)  Wx 

where  <Wt>  »  z  — ^ 

T-1  g(s) 


(24] 


(25] 


Upon  further  expansion  of  the  logarithm  above,  U(s)  can  be  written  as 


1  1 


U(s)  »  <W^> - {<Wt  >  -  <Wt>  } 

2!  KT 


+ - r  {<Wx^  -  3<Wt><Wt^  +  2<Wt>^},  (26] 

3!  (KT)^ 

where  <Wx^>  -  <Wt>^  *  <(Wx  -  <  W^>)^>  Is  the  second  moment  and  <W^^>  -  3  <Wx> 
<Wx^>  +  2  <Wt>^  -  <(Wt  -  <Wt>)^>  Is  the  third. 

Substituting  <Wt>  for  the  alloy  Fa  ■  Fg  ■  Into  the  above  equation  gives 
U(s), 
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r.-  ■.*  V  %■ 

■.  -•  .'v 


U(s)  -  —  Vab  +  V2N2V{1/2(1  -  s^)  - - (l-s^) 

2  KT  8 


V  2  1 


-  s"  (1-s^)^  +  — 

KT  12 


The  free  energy  Is  found  by  substituting  the  above  equation  and  equation 


(5)  for  g(s),  into  equation  (11), 


F(s)  -  F(l)  •  -  KT{(l+s)in(l+s)+(l-s)«n(l-s) 
2 


N  ,  V  1 


-  2tn2}  +  -  zV{V4(l-s  ) - (1-s  ) 

2  KT  8 


V  1 


- s^(l-s)^  + - }. 

KT  12 


The  B>W  solution  corresponds  to  keeping  only  the  first  and  second  terms  of  the 


above  equation  and  ignoring  a11  higher  order  terms.  Again  the  equilibrium 


value  of  s  at  T  is  found  by  minimizing  F(s), 


1+s  zV  VI 

An - »  —  s{l- - (1-5^) 

1-s  KT  KT  2 


V  2  1 


+ - (l-s^)(l-3s^)  +  — 

KT  6 


This  function  drops  off  much  more  rapidly  at  Tg  compared  to  the  B-M  approxima¬ 


tion  and  is  almost  identical  to  Bethe's  and  Guggenheims  quasichemical  model 


[26].  However,  difference  in  the  specific  heats  do  arise.  The  critical  temp- 


3F  3  F 


erature  is  again  found  by  the  condition  of  — 


s  ■  0,  and  is 


3S  3S 


KTc  IV  IV 

zV  2  KTc  6  KTc 


+ 


*  *  •_  I  •  •  ^  ^  ^  •_  •-  •-.m  m  m-  ^  - 


Kirkwoods  solution  has  the  advantages  over  Bragg  and  Williams  In  that  It  Is 
easily  adopted  to  different  types  of  lattices  and  to  nonstoichiometry.  Also, 
In  principle.  It  Is  possible  to  obtain  any  degree  of  accuracy;  however, 
equation  (30)  tends  to  converge  very  slowly  at  temperatures  below  and 
higher  moments  are  necessary  for  reliable  results. 

Bethe's  solution  [28]  has  been  shown  to  be  equivalent  to  Guggenheim's 
quasichemical  method  of  solution.  His  original  model,  as  will  be  presented 
here.  Is  basically  a  construction  of  the  grand  partition  function  for  a  small 
group  of  lattice  sites  In  which  the  variables  of  temperature,  volume  and  abso¬ 
lute  activities  play  Important  roles.  However,  the  original  model  was  not 
constructed  as  a  grand  canonical  ensemble.  Instead  Bethe  based  his  development 
on  nearest  neighbor  Interactions  and  the  short  range  order  parameter  a.  The 
short  range  order  parameter,  as  defined  by  Bethe,  Is  not  based  on  a  and  g 
sites  In  the  lattice  but  on  the  configuration  of  nearest  neighbors.  This 
parameter  measures  the  way  In  which,  on  the  average,  each  atom  Is  surrounded 
by  Its  neighbors,  or  the  extent  of  local  order.  The  short  range  order  Is 
defined  as 


(q  -  dr) 


(31) 


dm-dr 

where  q  »  Qab/Q>  dm  maximum  value  of  q  and  qp  Is  the  value  of  q  In  the 
disordered  state.  Here,  the  limits  of  zero  and  one  are  obtained  for  complete 
disorder  and  order  respectively,  as  for  the  long  range  order. 


Now,  we  select  an  arbitrary  group  of  sites  for  consideration  on  the  basis 


of  nearest  neighbor  interactions,  consider  for  example  Figure  2 


Figure  2;  Lattice  collection  in  Bethe's  approximations 

The  group  consists  of  an  interior  and  a  boundary,  all  remaining  sites  of  the 
lattice  are  exterior.  The  influence  of  the  exterior  atoms  are  taken  into 
account  only  on  the  average.  For  simplicity  the  method  is  developed  for  an  AB 
alloy. 

For  the  first  approximation  we  select  an  a  site  as  the  interior.  Thus  the 
boundary  will  consist  of  z  nearest  neighbor  s  sites.  For  the  second  approxi* 
mation  the  group  may  be  selected  as  in  Fig.  (2b);  the  larger  the  group  chosen, 
the  greater  the  accuracy,  however  the  mathematics  become  correspondingly  more 
complicated  [28].  Considering  the  first  approximations,  we  assume  that  a 
state  of  long  range  order  exists  in  the  exterior.  This  ordering  will  affect 
the  boundary  of  our  group  by  encouraging  A  atoms  to  a  sites  and  B  atoms  to  s 
sites.  This  tendency  to  order  on  the  boundary,  due  to  the  exterior  is  denoted 
by  an  ordering  energy  u.  The  energy  u  is  equal  to  the  difference  in  energy  of 
a  wrong  atom  on  an  a  or  3  site  and  a  right  atom  on  an  a  or  s  site. 


Consider  now,  the  energy  of  the  alloy  which  is  determined  by  the  number  of 

pairs  Qaa»  QbB  ®nd  Qab  for  a  definite  arrangement  of  atoms  on  the  lattice 

sites.  The  energy  of  a  crystal  of  pure  A  with  FaN  atoms  will  be  -  Vaa(#  of 
z 

pairs)  or  —  Fa^Vaa  *  "  FaQVaa  3od  similarly  for  a  crystal  of  B  atoms  only, 

2 

the  energy  is  -FbVbbQ*  Thus  the  energy  of  the  alloy,  with  respect  to  the 
pure  crystals  of  its  components  is 

E  ’  -  VaaQAA  -  VabQaB  -  VbbQbB  -  (“VaaFaQ  -  VbbFbQ)*  (32) 
This  reduces  to 

E  *  {VaA  +  Vbb)V!2  -  VabIQaB  *  -  VQab  (33) 

since 


(34) 


QAA  •  FaQ  -  Qab/2 
QBB  *  FbQ  -  Qab/2. 

So  that  the  superlattice  will  be  stable  at  lower  temperatures  V  must  be  taken 
less  than  zero.  The  Boltzmann  factor  is  then 

e-V/KT  ,  e  .  (35) 

For  the  boundary,  the  relative  probability  of  finding  an  A  atom  instead  of  a  B 
atom  due  to  the  exterior  is 


e-u/KT  .X.  (36) 

From  the  two  Boltzmann  factors  above  various  probabilities  in  the  interior  and 
in  the  boundary  may  now  be  found. 

If  the  central  atom  is  A,  a  right  atom,  the  relative  probability  of 
finding  n  wrong  A  atoms  in  the  boundary,  Pr(n),  is 

/z\ 


Pr(n) 


eOxO, 


(37) 


If  the  center  atom 'is  a  B  atom  the  probability  of  finding  n  wrong  atoms  on  the 
boundary  is 
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Pw  (n) 


ngZ-n, 


(38) 


since  each  of  the  z-n  right  B  atoms  in  the  boundary  has  an  interaction  energy 
V  with  the  wrong  interior  B  atom.  The  binomial  coefficient(n)is  the  number 
of  ways  to  arranging  n  wrong  A  atoms  in  the  z  boundary  sites.  The  total  rela¬ 
tive  probability  for  the  interior  atom  being  correct  is 


ri  =  I  Pr(n)  *  (1  +  ex)^  » 
n*0 


(39) 


and  for  being  wrong  is 


z 

W1  *  I  Pw(n)  *  (e  +  x)^.  (40) 

n*0 

The  normalized  probabilities  are  then 

ri  wi 

'*a  - - »  and  Wa  - -  .  (41) 

r^+wi  n+wi 

It  is  also  necessary  to  calculate  the  relative  probability  of  sites  in  the 
boundary  being  wrongly  occupied.  This  is  the  average  number  of  wrong  atoms  in 
the  boundary  divided  by  z.  The  relative  probability  of  finding  n  wrong  atoms 
in  the  boundary  is  Pr(n)  +  Pw(n),  the  average  number  of  wrong  atoms  in  the 
boundary  is  calculated  as  follows: 


z 

<n>  »  zwg  ■  z  n(Pr(n)  +  Pv<(n)) 

n»0  , 

Z  (42) 

Z  (Pr(n)  +  Pw(n)) 
n*0 


here,  wg  is  the  normalized  probability  of  having  a  wrong  atom  in  the  boundary 
6  sites,  therefore 


ex  X 

-  ri  + - 

1  +  ex  e+x 

wg  - -  .  (43) 

ri  +  Wi 


21 


Now  in  order  to  determine  the  unknown  energy  u  in  x  Bethe  introduced  a  sym¬ 
metry  or  consistency  condition.  Since  there  are  an  equivalent  number  of  o  and 
0  sites  and  A  and  B  atoms  for  an  AB  alloy  at  stoichiometry,  the  probability 
must  be  symmetric,  i.e.  the  same  probabilities  hold  whether  an  a  or  0  site  is 
chosen  as  the  origin  of  our  group.  Thus, 


Wa  “  (44) 
This  is  called  the  consistency  equation,  putting  equations  (41)  and  (43)  into 
the  above  equation  the  value  of  u  at  T  can  be  found  as  a  function  of  V  and  T, 


X  +  e 
1+xe 


(45) 


The  short  range  order  a  and  the  configurational  energy  of  the  alloy  can 
now  be  found.  The  number  of  AB  pairs  is  z-n  for  a  right  center  atom,  and  n 
for  a  wrong  center  atom,  therefore  on  the  average  the  fraction  of  AB  atoms  in 
the  alloy  is 


z  z 

z  (z  -  n)Pr(n)  +  i  nPw(n) 

QAB  n*0  n*0 

q  *  —  =  -  ^ 

Q  z(ri  +  wi ) 

For  the  AB  alloy  q^^^  »  1  and  q^  =  .5,  thus,  from  eq,  31 ,  cr  =  2q 
equation  (46), 

4ex  1 

a  »  1 - - — —  . 

1+ex  1+x^^'^”^' 

The  configurational  energy  is  then 

VQ 

E  »  VQaB  »  -  VQq  »  -  (ff  +  1), 
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measured  from  the  completely  ordered  alloy  this  energy  becomes 

VQ 

E  -  E(0)  - - (1-a). 


(46) 
or  from 


(47) 


(48) 


The  long  range  order  can  be  found  by  use  of  equations  (2),  (39),  (40),  and 
(41)  and  is 

(1+ex)^  -(e+x)^ 

s  *  2roi  -  1  - -  .  (50) 

(1+ex)^  +(e+x)^ 

Recalling  that  at  complete  disorder 

fo  =  fg  =  Wo  »  wg  »  V2 

and  using  equation  (45),  the  critical  temperatures  for  an  AB  alloy  Is  given  by 

- a  ___  {26]. 

KTc  z-2 

In  this  form,  extension  to  other  lattices  besides  simple  cubic  AB  alloys  Is 
quite  laborious.  However,  this  Is  equivalent  to  the  quasi  chemical  method  of 
Fowler  and  Guggenheim  [31];  and  Takagl's  combinatorial  analysis  [32]  allows  a 
simpler  extension  to  different  compositions  and  structure. 

A  comparison  of  the  three  early  models  can  be  seen  In  Figures  3,  4  and  5. 
All  three  of  these  models  have  neglected  Important  physical  considerations,  as 
do  the  more  recent  models,  these  Include:  (1)  Interaction  of  atoms  which  are 
not  nearest  neighbors,  (2)  changes  In  cell  size  and  lattice  symmetry,  and 
thermal  expansion,  (3)  lattice  vibrations,  (4)  limit  of  the  validity  of  the 
ordering  energy  and  finally,  (5)  effect  of  size  differential  In  the  ordering 
phenomena.  A  comparison  of  theory  and  experiment.  Figures  6  and  7  shows  poor 
quantitative  agreement,  qualitatively  the  models  are  essentially  correct  [26]. 
The  more  recent  models  are  somewhat  more  quantitatively  correct  and  will  be 
briefly  explained. 
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Figure  4:  Configurational  Energy  as  a  function  of  T.  [26] 


Figure  5:  Long-range  order  S  as  a  function  of  T.  [26] 


iii.  Cluster-Variation,  Concentration  Waves  and  Band  Theory 

The  cluster-variation  [39]  method  is  an  extension  of  the  quasi chemical 
model.  This  model  uses  the  same  type  of  energy  terms  as  the  earlier  theories; 
however,  it  introduces  a  way  of  arriving  at  a  better  approximation  to  the 
entropy.  The  free  energy  is  written  in  terms  of  the  different  configurations 
a  cluster  can  have  and  is  then  minimized  with  respect  to  the  configurations. 
For  example,  for  the  FCC  structure  a  tetragonal  group  of  atoms  is  chosen  as 
the  basic  cluster;  distribution  variables  are  assigned  to  each  configuration 
of  the  group  and  the  free  energy  is  written  in  terms  of  these  variables. 
Minimizing  the  free  energy  with  respect  to  these  variables  produces  a  set  of 
simultaneous  algebraic  equations,  usually  of  higher  order.  The  values  of  the 
distribution  variables  which  satisfy  these  equations  represent  the 
equilibrium  configurations  of  the  system  [40].  The  early  formulation  of  the 
method  did  not  allow  one  to  choose  very  large  clusters  (which  increases  the 
accuracy)  due  to  the  large  number  and  nature  of  the  simultaneous  equations. 
Kikuchi,  however,  has  developed  an  iteration  scheme.  Natural  Iteration  [39], 
which  is  particularly  suited  to  the  method  and  guarantees  convergence. 

A  number  of  authors,  [41,  42,  43],  have  used  this  method  and  results 
agree  fairly  well  with  experiment  and  have  been  reproduced.  This  is  perhaps 
one  of  the  more  promising  methods  of  describing  cooperative  phenomena. 

The  idea  of  concentration  waves  was  introduced  by  Landau  [44]  and  Lifshitz 
[45]  in  (1937,  1941)  in  the  development  of  the  phenomenological  theory  of 
phase  transitions  of  the  second  kind.  This  method,  applied  to  the  order- 
disorder  transition  in  alloys  [46]  has  proved  to  be  fruitful.  This  theory 
enables  one  to  understand  the  symmetry  aspects  of  the  order-disorder  tran¬ 
sitions,  to  take  into  consideration  long  range  interactions  in  the  alloy  and 
to  predict  the  structure  of  the  ordered  phases,  which  the  other  theories  do 

not  allow,  based  on  concrete  interatomic  potentials. 
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All  of  the  preceedlng  theories  reduce  to  the  comparison  of  the  free  energy 
of  the  ordered  and  disordered  phases  calculated  from  nearest  or  at  best  next 
nearest  Interactions.  The  question  of  whether  or  not  the  structure  of  the 
ordered  phase  Is  stable  [at  the  relevant  choice  of  Interatomic  energies  I.e. 
what  value  of  1  Is  chosen  In  equation  (1)]  and  what  the  structure  of  the 
ordered  phase  Is  are  not  considered.  The  early  theories  as  we11  as  the 
cluster-variation  Improvement  become  especially  unsatisfactory  In  light  of 
recent  developments  In  the  pseudopotential  theory  of  metal  alloys  and  cova¬ 
lent  compounds  (a  large  part  of  the  T1-C  bonds  are  covalent)  as  well  as  the 
theory  of  stress  Induced  Interactions  In  a  solid  solution.  Studies  show  that 
Interatomic  Interactions  which  are  responsible  for  ordering  cannot  be  short 
range.  Thus  for  substantial  Improvement  of  the  above  theories  the  real  long 
range  Interactons  must  be  taken  Into  account. 

These  difficulties  are  overcome  by  using  the  method  of  static  con¬ 
centration  waves.  Atomic  distribution  In  a  binary  alloy  can  be  described  by 
means  of  the  function  n(7n*  this  one  function  Is  adequate  because  the  occupa¬ 
tional  probabilities  nA("r)  and  nBCr)  for  A  and  B  atoms  respectively  are  not 
Independent,  they  must  satisfy  the  Identify 

haCH  +  neCr)  =1,  (52) 

here  7  Is  a  crystal  lattice  vector.  This  function  can  be  used  for  Intersti¬ 
tial  as  well  as  substitutional  alloys.  In  the  disordered  state  the  probabili¬ 
ties  n(7)  are  the  same  for  all  sites  which  can  be  occupied.  They  are  equal  to 
the  atomic  fraction  c  of  the  relevant  component  In  a  substitutional  alloy  and 
In  an  Interstitial  solution  they  equal  the  fraction  of  Interstitial  sites 
which  are  occupied.  In  the  ordered  phases  n(7)  becomes  dependent  on  the  site 
coordinate  7.  For  Instance,  If  n(7)  describes  t  values  ni,  n2 — n^  on  a  set 


of  crystal  sites  {7},  then  nCr)  describes  the  sublattices  Into  which  the 
disordered  lattice  has  been  subdivided  as  a  result  of  ordering.  The  n^  are 
occupational  probabilities  of  the  sites  of  the  1st,  2nd  —  t-th  sublattice. 

The  function  n(7)  can  be  expanded  In  the  form  of  a  Fourier  series,  I.e.  It  can 
be  represented  as  a  superposition  of  static  concentration  waves: 

n(r)  «  C+V2  ?CQ(kj)e^*^J^  +  (}*(Trj)e“'*l^j^»  ^53) 

where  exp(lkr)  Is  a  static  concentration  wave,  Tj  Is  a  non-zero  wave  vector 
defined  In  the  first  Brillouin  zone  of  the  disordred  alloy,  7  Is  a  site  vec¬ 
tor  of  the  lattice  {7}  describing  a  position  which  can  be  occupied  by  an 
atom,  the  Index  j  denotes  the  wave  vectors  In  the  Brillouin  zone  and  QCEj)  Is 
a  static  concentration  wave  amplitude.  Alternatively  this  can  be  written  by 
combining  those  terms  whose  wave  vectors  TTj's  enter  into  the  stars  (the  star 
Is  a  set  of  wave  vectors  ITj,  which  may  be  obtained  from  one  wave  vector  by 
applying  to  It  all  operations  of  the  symmetry  group  of  the  disordered  phase): 

n(7)  »  c  +  E  n^E  (7),  (54) 

s  * 

where  Es(7)  -  1^  s(v(js)exp(fEjj7)  +  (55) 

Js 

and  QCEjj)  «  nsTrs{js)»  (56) 

The  summation  Is  carried  out  over  all  wave  vectors  In  the  star,  ns  are  the 
long  range  order  parameters  and  ys(Js)  are  coefficients  which  determine  the 
symmetry  of  the  occupation  probabilities  n(7)  with  respect  to  rotation  and 
reflection  symmetry  conditions.  The  requirements  that  In  the  completely 
ordered  state,  when  all  n(7)  are  either  zero  or  unity,  that  all  nj  are  equal 
to  unity  completely  defines  the  constants  Ys(js)*  It  turns  out  that  the  con¬ 
centration  waves  are  related  to  the  amplitude  of  coherent  scattering  In  X-Ray 
diffraction  by  the  equation 


(57) 


Y(q)  -  (fB-fA)V!2i:CQ(kj)5:e-i(q-kj)r 
+  Q*(kj)ze-l(q+kj)r3  , 

r 

where  and  fg  are  the  form  factors  of  the  A  and  B  atoms  respectively. 

Now*  the  configuration  part  of  the  Hamiltonian  of  the  alloy  is 

H  »  V!2  £  V(r,7")c(r)cCr'),  (58) 

rr' 

where  the  summation  is  over  all  of  the  lattice  sites  and 

c(F)  *  1  if  there  is  a  solute  atom  at  r, 
c(7)  *  0  if  not. 

The  interchange  energy  V(7,F')  is  defined  by  the  equation 

V(r.r")  »  VAA(r.r")  +  VagCr.?')  -  2VAB(r.r').  (59) 

This  is  analogous  to  the  energy  defined  in  equation  (1),  however  it  takes  into 
account  all  interactions  in  the  lattice. 

The  problem  of  knowing  the  ordered  phase  structure  reduces  to  determining 
the  function  n(?).  Notice  that  n(r)  follows  a  sort  of  Pauli  exclusion  prin¬ 
ciple:  each  lattice  site  is  either  occupied  by  one  or  zero  atoms  of  some 

definite  type.  Thus  we  can  reason  that  it  follows  the  Fermi -Dirac 
distribution: 

1 

nCr)  - -  ,  (60) 

-p+*(7) 

exp - +1 

KT 


where  the  chemical  potential  u  is  determined  by  the  condition  of  conservation 
of  mass  t—  ^ 


1 


zn(r)  ■  E 
r  r 


-u+<fr(r) 

KT 


•  Ni 


(61) 
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(62) 


where  Is  the  number  of  solute  atoms.  In  the  self-consistent  field  approxi 
mation,  ^{7)  is 

♦  (7)  -  E  V{7,?')n(?). 

This  approximation  neglects  correlation  effects. 

Substituting  this  potential  into  (60) 

1 


n(r) 


-u  1  _ _ 

exp  —  —  E  V(r,r')n(r')  +1 

KT  KT  r 


(63) 


which  corresponds  to  the  free  energy 


n  «  U-TS  -  u2n(7)  , 
r 


(64) 


where  U  ■  1/2  £  VCr,r')n(7)n(r‘’) 
r,r' 


is  the  internal  energy  and 


S  »-kE[n(r)An(n(T^)  -*•  (1  -  n(7))4n(l-n(7))] 
r 


(65) 


is  the  entropy.  The  Helmholtz  free  energy  of  the  system  is 

F  -  U-TS. 

By  using  the  expression  for  n(7)  in  equation  (54),  and  equation  (63)  the 
equilibrium  values  of  ns  at  T  can  be  found: 


(66) 


t-1 


c  +  E  nsEs(^  *  ®xp 
s 


-w+v(0)c+  E  V(rs)nsEs('fO 
s-1 


+  1. 


KT 


(67) 


Here  VCIT)  «  EVCr)  is  the  Fourier  transform  of  the  interatomic  potential, 
and  can  be  found  by  measuring  the  diffuse  scattering  intensity  from  a  single 
crystal  of  the  disordered  alloy,  [46] 

c(l-c) 

(68) 


IdiffC^  ■  (^A  -  ^B)^ 


c(l-c) 

1  + - v(ir) 


kT 
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This  method  can  predict  the  critical  temperature,  ordered  structure  and  stabi¬ 
lity  apparently  with  much  more  ease  and  accuracy  than  the  former  theories 
[46].  The  critical  temperature  calculated  for  the  CuAu3  alloy  by  this  method 
Is  only  19**  lower  than  the  actual  experimental  value.  This  accuracy  Is 
exceptional  as  compared  to  the  previous  models;  It  also  has  the  advantage 

I 

that  the  Interaction  energies  are  In  no  way  ambiguous  and  are  measurable 
through  X-Ray  diffraction. 

One  of  the  band  theory  models  of  the  ordering  phenomena  [37]  uses  a  tight 
binding  Hartree  Hamiltonian  to  determine  the  ordering  energy  Incorporated 
with  the  cluster-variation  method  for  the  entropy  of  the  system.  Here  the 
bulk  density  of  states  Is  derived  from  a  configurational  average  over  the 
local  densities  of  states  at  the  central  site  of  a  cluster.  The  free  energy 
of  the  alloy  Is  determined  by  minimizing  a  model  free-energy  function  over  the 
space  of  possible  electron-ion  configurations.  These  configurations  may  be 
classified  by  the  values  of  certain  order  parameters  which  measure  multi - 
site  correlations,  e.g.,  long-range  order,  short  range  order,  three-body 


correlations,  etc.  Thus  at  a  temperature  T, 

F(T,c)  »  mln^j^  F(T,c;oil) 

=  E{T,c;aj^)-TS{c,aj^) 


(69) 


where  c  Is  the  composition.  A  schematic  of  the  free  energy  calculation  Is 


shown  In  Figure  8. 
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In  Figure  8.,  Nina  M  the  local  density  of  states  and  <nino>  is  the  ther¬ 
mal  averaged  occupation  of  sites. 

This  model  predicts  a  critical  temperature  of  the  CuAu  system  which  is  in 
error  by  19%  which  is  more  accurate  than  the  early  set  of  models  but  can  not 
compare  with  the  accuracy  obtained  in  the  method  of  concentration  waves.  One 
notable  point  of  this  theory  however  is  that  it  actually  calculates  the  atomic 
interaction  energies  as  a  function  of  short-range  order,  these  energies  are 
constant  only  near  a*0  and  can  vary  in  magnitude  as  much  as  60%.  This  makes 
the  early  Ising  models  and  the  cluster-variation  model  even  more  unsatisfac¬ 
tory  since  they  are  based  on  an  a  priori  knowledge  of  this  energy,  which  is 
assumed  constant  at  all  values  of  order. 

Since  Fosdick  [38]  first  applied  the  Monte  Carlo  Method  to  order-disorder 
transformations  in  solids  a  number  of  attempts  have  been  made  to  improve  and 
extend  this  procedure  [47].  Unfortunately  the  stability  of  any  given  con¬ 
figuration  in  this  model  is  again  based  on  the  nearest  neighbor  interaction 
energies.  Since  the  method  can  require  quite  a  lot  of  computation  time  and 
exhibits  the  same  energetic  restrictions  as  early  Ising  theories  this  method 
of  approaching  the  problem  of  ordering  does  not  seem  to  be  worth  the  efforts 
in  our  present  research  effort. 


1v.  Ordering  Energy 

It  seems  that  the  most  productive  model  available  to  the  predictions  of 
the  ordering  phenomena  lies  in  the  method  of  static  concentration  waves  and 
this  is  probably  the  theory  which  we  will  use  in  predicting  the  ordering  pro¬ 
perties  of  the  Ti-C  system.  The  absence  of  ambiguity  as  far  as  energy  terms 
are  concerned  makes  this  theory  much  more  desirable  than  the  Ising  models. 

The  band  model  is  orders  of  magnitude  more  involved,  as  well  as  less  accurate. 
However,  it  will  be  desirable  to  compare  predictions  made  by  the  Ising  models 
and  to  this  end  an  understanding  of  the  ordering  energy  in  equation  (1)  is 
necessary. 

Some  Justification  of  choosing  nearest  neighbor  interaction  energies 
Instead  of  long-range  effects  in  ordering  seems  to  be  found  in  the  quantum 
theory  which  holds  that  forces  between  atoms  decrease  rapidly  as  the  distance 
separating  them  increases.  As  stated  earlier,  the  pair  Interaction  energies 
Vaa>  interaction  energies  between  pairs  of  atoms  in  the  pure  ele¬ 

mental  crystals  and  Vab  the  interaction  in  the  alloy.  However,  this  defini¬ 
tion  is  very  ambiguous.  What  do  we  measure  to  find  these  energies,  or  how  do 
we  calculate  them?  The  energy  in  equation  (1)  is  the  first  term  in  the 
expansion  of  the  Fourier  Transform  of  the  long-range  interaction  energy  used 
in  the  static  concentration  wave  model,  where  for  a  FCC  lattice,  [46] 

V(Ko)  “  V(2Tra*3)  »  -4wi  +  6w2  -  8W3-  —  (70) 

where  w^  are  the  energies  corresponding  to  each  coordination  shell. 

The  simple  nearest  neighbor  approximation  has  been  shown  to  be  extremely 
crude  [49,50],  for  instance,  it  cannot  explain  why  most  order-disorder 
transformations  are  first  order,  why  the  maximum  value  of  T^  often  occurs 
away  from  stoichiometry  or  why  an  FCC  alloy  should  order  at  all.  Freidel  [50] 
and  others  [51]  have  been  able  to  show  that  because  conduction  electrons  in  a 
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metal  are  limited  to  energies  below  the  Fermi  surface  a  complete  screening  of 
the  icoic  metal  cores  are  not  possible,  thus  long  range  effects  are  present. 
The  pseudo-potential  theory  [52]  also  points  to  many  body  interactions. 
However,  neither  one  of  these  methods  of  calculating  interaction  energies  can 
be  used  for  the  alloy  system  under  investigation  with  any  faith  due  to  the 
presence  of  covalent  bonding.  These  bonds  have  a  strong  angular  nature  and 
many  body  interactions  [23]. 

A  theory  recently  developed  to  find  the  ordering  energies  of  transition 
metal  alloys  [53]  is  available,  however  it  has  not  been  used  and  requires 
quite  a  few  parameters  which  might  be  as  difficult  as  the  ordering  energies 
themselves  to  obtain. 

In  addition  to  the  extensive  literature  search  for  methods  to  calculate 
and  measure  the  ordering  energies  we  also  consulted  several  members  of  the 
faculty  in  chemistry  and  solid  state  physics  here  at  Auburn  about  this  problem 
and  without  exception  they  all  suggested  using  heats  of  formations  for  these 
energies.  After  looking  in  several  references  [54,55]  we  found  values  that 
were  experimentally  verified  (see  table  1)  and  used  these  in  calculating  the 
critical  temperature  for  the  Cu-Au  alloy  system. 


bond 

energy 

Au-Au 

88  ^0.5  kcal/mole 

Cu-Cu 

80.7  kcal/mole 

Au-Cu 

55.4  ^  2.2  kcal/mole 

Table  1.  Heats  of  formation  of  standard  states  from  monatomic  gases. 

The  critical  temperature  calculated  from  these  values  using  the  Bragg-Williams 
approximation  was  five  orders  of  magnitude  higher  than  the  measured  value. 
Clearly  the  heats  of  formation  are  not  the  interaction  energies  which  we  are 
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looking  for,  although  they  must  be  connected  In  some  way.  The  literature  had 
nothing  to  say  on  this  connection  however  and  the  idea  of  using  the  heats  of 
formation  for  ordering  energies  has  been  abandoned.  It  now  appears  that  the 
only  way  to  arrive  at  an  acceptable  number  for  these  energies  is  through  X-Ray 
diffraction. 


V.  Sunmary 

Although  there  have  been  quite  a  number  of  theories  put  forth  since  the 
pioneer  work  of  Bragg-Williams.  most  of  them  seem  to  exhibit  the  same  short¬ 
comings  (though  perhaps  not  as  severe)  as  this  first  approximation.  It  seems 
however,  that  the  use  of  the  method  of  static  concentration  waves,  coupled 
with  X-Ray  diffraction  experiments  will  give  satisfactory  results.  It  has 
been  somewhat  distressing  to  find  after  a  very  extensive  literature  search, 
(both  with  a  computer  and  in  the  library)  that  very  little  work  has  been  done 
in  actually  finding,  theoretically  or  experimentally,  the  needed  interaction 
energies.  Although  there  is  one  theory  available  to  calculate  the  ordering 
energies  [53],  no  one  has  used  it  and  it  is  doubtful  that  we  will  be  able  to 

due  to  its  complexity.  Clapp  and  Moss  have  done  some  diffraction  work  in  the 

Cu-Au  system  in  which  they  have  tried  to  find  the  pair-wise  interaction 
energies,  however,  basically  what  they  have  done  is  fit  this  to  the  known  cri¬ 
tical  temperatures  [57].  Semenouskaya  [58]  using  X-Ray  work  done  by  Moss  [59] 

has  attained  remarkable  results  in  calculating  the  ordering  temperature  of 
Cu3Au  as  stated  before  and  it  is  hoped  that  we  will  be  able  to  employ  this 
same  method  and  attain  the  same  kind  of  results.  The  measurements  made  were 
the  absolute  intensity  of  X-Ray  diffuse  scattering  by  the  disordered  alloy. 
There  are  a  few  references  in  the  literature  on  diffraction  studies  done  with 
single  crystal  Ti-C,  however  we  will  undoubtedly  require  more  data  than  is 
available.  To  this  end  we  will  have  to  take  diffraction  patterns  with 
powdered  samples. 

Now  that  the  literature  search  and  study  is  nearing  completion  the  main 
problem  to  be  tackled  at  this  time  is  the  calculation  or  measurement  of 
ordering  energies  so  that  predictions  of  the  critical  temperatures,  etc.  can 


be  made.  To  this  end  powder  samples  from  the  alloys  already  melted  will  be 
made  and  an  analysis  of  the  diffuse  scattering  of  X-Rays  from  the  samples 
carried  out.  From  the  Intensities.  It  Is  hoped  that  an  accurate  measurement 
of  the  ordering  energies  will  be  possible.  The  effort  to  theoretically  come 
by  these  energies  will  also  continue.  Once  these  energies  are  found  the 
ordering  properties  of  the  alloy  system  will  be  studied  at  various  com¬ 
positions  and  the  most  promising  of  these  compositions  will  be  concentrated  on 
In  the  experimental  phase. 

Besides  predicting  the  structure  and  ordering  temperature  of  the  alloys, 
compression  tests  at  elevated  temperatures  are  planned,  as  well  as  measure¬ 
ments  of  the  specific  heats.  Once  the  properties  of  these  alloys  are  well 
understood,  additional  elements  may  be  alloyed  with  the  T1-C  to  Improve  any 
undesirable  properties  such  as  brittleness. 

Efforts  to  determine  the  effects  of  additional  alloying  elements  will  be 
made  after  Kikuchl  [24].  These  calculations  In  conjunction  with  electron  den¬ 
sity  calculations  will  point  to  the  best  choice  of  alloy  additives.  These 
studies  will  form  a  strong  basis  for  the  additional  study  and  optimization  of 
the  alloy  system. 


B.  Experimental  Progress 

The  C-T1  system  was  Identified  as  the  most  promising  alloy  system  based 
upon  the  ability  to  obtain  an  ordered  structure  over  a  range  of  compositions 
and  the  ultrahigh  melting  temperature  of  the  system  (melting  point  of  3080  C). 
A  series  of  C-Ti  alloys  was  melted  to  investigate  the  properties  of  the 
system.  Fifteen  compositions  were  melted  covering  the  range  0  to  22  weight 
percent  carbon.  The  Initial  alloy  melts  did  not  include  additions  for  impro¬ 
vement  of  ductility  to  enable  a  direct  comparison  between  theoretical  predic¬ 
tions  of  ordering  parameters  and  experiment.  Initial  melts  were  made  by  arc 
melting  under  an  Argon  cover  gas. 

The  arc  melting  unit  used  to  produce  the  initial  C-Ti  melts  consists  of  a 
gas  cooled  electrode  (0.38  inch  diameter  stinger)  which  is  composed  of 
Tungsten-2%  Thoria  and  a  graphite  crucible  approximately  .5  inches  in  diameter 
by  .25  inches  high.  To  produce  these  melts,  direct  current  voltages  of  35-40 
volts  and  450-575  amps  were  required.  Some  problems  were  encountered  during 
the  arc  melting  process  because  of  the  high  melting  temperatures  needed  to 
produce  the  alloys.  An  inert  gas  flow  rate  that  was  higher  than  normal  was 
necessary  to  keep  the  Tungsten  stinger  from  melting.  This  inert  gas  flow  had 
a  tendency  to  blow  the  carbon  additions  out  of  the  crucible,  resulting  in  some 
of  the  melts  having  a  carbon  composition  less  than  desired.  Compositions  of 
the  alloys  have  not  been  checked  and  hence  reported  compositions  are  nominal 
compositions  which  do  not  reflect  such  losses.  Samples  produced  by  this 
method  were  restricted  to  buttons  less  than  5  cubic  centimeters  in  volume. 

For  melts  larger  than  this  it  was  impossible  to  keep  the  entire  specimen 
molten  throughout  the  procedure. 
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Figures  9  through  12  show  representative  microstructures  of  the  fifteen 
alloys  produced  by  this  arc  melting  procedure.  Figure  9  Is  the  photo- 
micrograph  of  a  1  weight  per  cent  (w/o)  carbon-titanium  alloy  which  was 
generated  by  the  arc  melting  procedure  In  a  carbon  crucible.  Figures  10 
through  12  are  microstructures  of  alloys  with  Increasing  carbon  content.  The 
figures  show  that  the  microstructures  that  one  obtains  become  more  complex  as 
carbon  content  Is  Increased.  In  the  low  carbon  alloy  (Figure  9)  the 
microstructure  Is  dominated  by  a  single  phase  with  carbide  particles  distri¬ 
buted  along  the  grain  boundaries.  However  the  22  weight  percent  carbon  alloy 
structure  (Figure  12)  consists  of  at  least  four  distinct  phases.  Different 
structure  morphologies  were  also  Identified  from  the  optical  micrographs.  The 
microstructure  of  Figure  10  (a  12  w/o  carbon  alloy)  is  dominated  by  a  dendritic 
structure  while  the  microstructure  of  Figure  11  (a  16  w/o  carbon  alloy  is 
dominated  by  a  Widmannstatten  pattern. 

Figure  13  Is  the  phase  diagram  for  the  C-T1  system.  Comparing  the  results 
with  published  phase  diagrams  for  the  C-T1  system  one  finds  good  agreement 
between  the  anticipated  phases  and  the  microstructure  for  alloys  of  12  w/o 
carbon  and  below.  It  should  be  pointed  out  that  the  phase  diagram  of  this 
system  has  not  been  extensively  Investigated.  For  carbon  concentrations 
greater  than  12  w/o  and  temperatures  above  1500**C,  the  diagram  Is  based  on 
theoretical  calculations.  For  compositions  with  carbon  content  less  than  12 
w/o,  the  structure  should  consist  of  a  mixture  of  a-T1  and  y  phases.  This  Is 
what  was  obtained  In  the  alloys  shown  In  Figure  9  (1  w/o  carbon)  and  Figure  10 
(12  w/o  carbon).  However  for  alloys  with  carbon  concentrations  above  12  w/o, 
more  than  two  phases  were  found  In  the  microstructures.  According  to  the 
phase  diagram,  the  phases  present  should  be  a  or  Y+oraphlte.  One  of  the  pha¬ 
ses  In  Figure  11  and  Figure  12  Is  graphite,  but  there  appears  to  be  at  least 
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Ti  plus  1  w/o  Carbon  Alloy,  Figure  10:  Ti  plus  12  w/o  Carbon  Alloy. 

Structure  consists  of  a-Ti  phase  The  alloy  structure  consists  of  y 

with  TiC  carbides  at  grain  boundaries.  matrix  with  o-Ti  phase  dendrites. 


two  other  phases  which  do  not  appear  on  the  phase  diagram.  X-ray  diffraction 
studies  are  being  performed  In  an  attempt  to  Identify  the  phases  present  In 
the  alloys  of  Figure  10  and  12,  as  well  as  the  degree  of  order  achieved  In  the 
y  phase. 

Hardness  measurements  were  made  on  the  alloys  to  scope  the  mechanical  pro¬ 
perties.  The  macroscopic  hardness  measurements  were  found  to  be  relatively 
Insensitive  to  carbon  content.  The  Rockwell  C  hardness  varied  from  a  value  of 
20  for  the  1  w/o  carbon  alloy  to  52  for  the  22  w/o  carbon  alloy.  The  high 
carbon  alloys  were  found  to  crack  If  the  hardness  Indentation  was  applied  near 
the  edge  of  the  test  specimen.  In  addition  to  macrohardness  measurements,  an 
attempt  to  make  microhardness  measurements  on  the  various  phase  structures 
seen  In  Figures  9  through  12  was  made.  This  attempt  was  only  partially  suc¬ 
cessful.  Only  on  the  o-T1  phase  could  valid  readings  be  obtained.  Attempts 
to  make  Indentations  on  other  phase  structures  were  unsuccessful .  No  detec¬ 
table  Indentations  occurred  for  the  applied  loads.  It  has  been  hypothesized 
that  the  other  phases  have  a  hardness  greater  than  DPH  2000  which  Is  the  upper 
limit  of  our  microhardness  test  machine  at  the  present  time.  A  special  weight 
pan  Is  being  purchased  to  Increase  the  load  that  can  be  applied  to  the  speci¬ 
men  and  microhardness  tests  will  again  be  made. 

Recently  a  new  Induction  heating  system  was  delivered  to  Auburn  which  was 
purchased  under  this  contract.  It  consists  of  a  5  kilowatt  generator  with 
Induction  coll  to  which  a  quartz  tube  for  maintaining  an  Inert  gas  atmosphere 
has  been  attached.  About  a  month  of  effort  was  spent  tuning  the  frequency  of 
the  Induction  generator  and  experimenting  with  coll  design  to  obtain  the 
higher  temperatures  required  In  melting  and  sintering  the  carbon-titanium 
alloys.  The  coll  shape  and  diameter  of  the  spiral  controls  the  concentration 
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Figure  13:  Phase  diagram  for  the  C-TI  alloy  system. 
Note  that  the  dashed  lines  indicate 
calculated  phase  relations. 


of  energy  on  the  melting  volume.  The  frequency  of  the  Induction  generator  can 
be  varied  to  Induce  deep  penetration  or  surface  heating  of  the  material  or 
maximize  the  heating  effect  as  a  function  of  alloy  system.  Although  higher 
melting  temperatures  are  obtained  as  the  diameter  of  coll  spiral  Is  decreased, 
the  size  of  the  melt  and  the  fabrication  of  the  coll  becomes  more  difficult. 
The  diameter  of  the  copper  tube  Is  fixed  In  the  present  system  In  order  to 
maintain  adequate  cooling  of  the  Induction  coll.  A  photograph  of  this  system 
In  operation  Is  shown  In  figure  14. 

With  this  simple  system  an  additional  ten  melts  have  been  made  of  carbon- 
titanium  alloys  covering  the  range  0-14  weight  per  cent  carbon.  Although 
there  are  still  problems  with  obtaining  a  homogeneous  material  composition  and 
eliminating  porosity  of  the  melt,  alloys  produced  using  this  system  are 
greatly  Improved  over  those  produced  by  previous  methods.  We  continue  to 
experience  difficulty  In  obtaining  a  fluid  mixture  of  the  alloy  for  carbon 
compositions  above  10  w/o.  Sintering  of  the  carbon  and  titanium  Is  obtained 
but  there  Is  no  real  melting.  From  these  newly  produced  alloys  two 
compression  samples  were  cut  for  mechanical  properties  testing.  The  carbon 
content  of  these  titanium  alloys  were  4  w/o  and  5  w/o  carbon.  Macrohardness 
tests  were  conducted  on  the  alloys  yielding  a  value  between  30  and  40  Rockwell 
C. 

A  diamond  slitting  saw  was  used  to  prepare  cube  compression  specimens 
because  of  the  small  size  of  the  available  material.  Although  the  alloys  were 
not  excessively  hard,  large  amounts  of  time  were  required  to  prepare  the 
compression  samples  which  measured  2.29mm  X  2.29mm  X  12.7mm  high.  The  alloys 
were  gummy  In  nature  requiring  sharpening  of  the  diamond  blade  with  a  carbon 
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block  to  remove  cut  particles  stuck  In  the  blade.  Compression  tests  were  con¬ 
ducted  at  room  temperature  using  a  Materials  Testing  System.  The  data  were 
simultaneously  recorded  on  a  HP  X-Y  chart  recorder  and  through  a  digital  HP 
9836  data  acquisition  system.  Figure  15  Is  the  stress-strain  curve  for  the 
compression  test  of  the  5  w/o  carbon  sample. 

A  summary  of  the  mechanical  properties  of  the  two  titanium-carbon  alloys 
along  with  properties  for  pure  titanium  and  T1-6A1-4V  alloys  Is  given  1n  Table 
2.  The  pure  T1  and  T1-6A1-4V  data  was  taken  from  the  ASM  Metals  Handbook  [60] 
and  Is  based  on  tensile  specimens.  Both  of  the  compression  samples  showed 
visible  signs  of  porosity  on  the  surface  prior  to  testing  (Figure  16). 

Despite  the  observable  flaws  In  the  specimens,  exceptional  yield,  ultimate  and 
elongation  values  were  obtained  for  the  5  w/o  carbon  alloy.  A  yield  strength 
of  100  ksl,  UTS  of  164  ksl  and  elongation  of  5.4%  Is  very  Impressive  when  com¬ 
pared  to  a  120  ksl  yield  strength,  UTS  of  130  ksl  and  4%  elongation  of  forged 
T1-6A1-4V.  In  addition  the  melting  temperature  of  the  5  w/o  carbon  alloy  Is 
expected  be  2000®C  compared  to  the  melting  point  of  T1-6A1-4V  of  1650®C. 

Both  compression  samples  failed  on  a  plane  of  approximately  45  degrees  to 
the  tensile  axis  along  a  line  connecting  several  of  the  pores  on  the  speci¬ 
men's  surface.  The  4  w/o  carbon  specimen  failed  prematurely  due  to  an 
Interior  gas  pocket  which  was  visible  after  the  test.  Optical  photo¬ 
micrographs  of  the  two  structures  are  given  In  figures  17  and  18.  Two 
separate  distinct  phases  appear  which  Is  consistent  with  the  phase  diagram. 

The  properties  of  the  alloys  appear  to  be  reversed,  however,  the  5  w/o  C  alloy 
Is  less  hard  than  the  4  w/o  alloy.  Verification  of  the  alloy  composition 
needs  to  be  confirmed  prior  to  any  conclusions  being  drawn  as  to  the  effect  of 
Increasing  carbon  content  on  the  mechanical  properties. 


Table  I 


Summary  of  Mechanical  Properties 


Alloy 

Melting 

Temp(®C) 

Yield 

(Ksi) 

UTS 

(Ksi) 

Ti-4C 

1900 

42.4 

54.2 

Ti-5C 

2000 

100 

164 

Ti 

1665 

20.3 

34.1 

Ti-6A1-4V 

1650 

120 

130 

El ongation 
(Ksi) 

Modulus 

(loVi) 

Hardness 

1.4% 

12.8 

41Rc 

5.4% 

30.8 

32Rc 

54% 

16.8 

70  HB 

4.0% 

15.4 

32Rc 

Fig.  14  Induction  generator  and  melting  furnace  used  to  produce  Ti-C  alloys 
for  compression  testing 
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Based  upon  our  Initial  experience  with  the  melting  of  the  C-Ti  alloys,  two 
of  the  principle  areas  which  must  be  improved  are  grain  boundary  structure  and 
uniformity  of  composition  throughout  the  melt.  In  alloys  with  carbon  addi¬ 
tions  greater  than  12  w/o,  the  graphite  is  found  to  reside  along  the  grain 
boundaries  or  as  blocky  phases  in  the  matrix.  A  more  uniform  distribution  of 
constituents  is  anticipated  to  improve  the  ductility  and  also  enhance  the 
strength.  An  additional  improvement  which  is  being  pursued  is  the  elimination 
of  grain  boundary  impurity  elements  such  as  the  halide  elements  which  are 
naturally  attracted  to  the  grain  boundaries.  An  ultrahigh  vacuum  furnace 
capable  of  2700**C  is  being  assembled  for  melting  and  heat  treatment  of  the 
second  generation  of  C-Ti  alloys.  The  future  melts  of  C-Ti  alloys  will  con¬ 
tain  between  1  and  5  w/o  Nb  and  0.5  w/o  Ce  to  improve  the  matrix  ductility 
through  control  of  the  e/a  ratio  and  reduce  the  segregation  of  impurity  ele¬ 
ments  at  the  grain  boundaries. 


V.  CONCLUSIONS 

•  A  theoretical  study  to  predict  critical  ordering  temperatures  and  com¬ 
positions  In  ordered  alloys  has  been  initiated.  Scoping  calculations  have 
Identified  the  C-T1  system  as  the  most  promising  system  In  which  the  ultra- 
high  strength,  ductility  and  temperature  resistance  can  be  obtained. 

•  A  total  of  twenty-five  T1-C  alloys  have  been  manufactured  covering  com¬ 
positions  of  0-22  w/o  C.  Optical  metallography  on  specimens  with  greater 
than  12  w/o  seem  to  show  phases  which  do  not  appear  on  standard  phase 
diagrams  of  the  C-T1  system. 

•  Compression  tests  on  a  5  w/o  C  sample  at  room  temperature  show  a  yield 
stress  of  100  ksl,  ultimate  tensile  strength  of  164  ksl  and  strain  to  frac¬ 
ture  of  5.4%.  These  properties  are  better  than  forged  bar  T1-6A1-4V 
properties. 
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